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Abstract 

We investigate electron shuttling in three-terminal nanoelectromechanocal device built on a mov¬ 
able metallic rod oscillating between two drains. The device shows a double-well shaped electrome¬ 
chanical potential tunable by a source-drain bias voltage. Four stationary regimes controllable by 
the bias are found for this device: (i) single stable fixed point, (ii) two stable fixed points, (iii) 
two limiting cycles, and (iv) single limiting cycle. In the presence of perpendicular magnetic field 
the Lorentz force makes possible switching from one electromechanical state to another. The 
mechanism of tunable transitions between various stable regimes based on the interplay between 
voltage controlled electromechanical instability and magnetically controlled switching is suggested. 
The switching phenomenon is implemented for achieving both a reliable active current switch and 
sensoring of small variations of magnetic field. 

PACS numbers: 81.07.Oj, 73.23.Hk 
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Nanoelectromechanical (NEM) systems arouse interest not only due to diverse potential 
applications as nano-devices but also as an efficient constituent of modern nano-electronics 
DEI- While the NEM coupling plays an important part in electronic transport through 
nano-devices, the charge transport associated with the nano-mechanical motion demon¬ 
strates various interesting quantum effects, such as Coulomb blockade |3], resonant tunnel¬ 
ing @J, spin-dependent transport (5] and so on. Besides, strong NEM coupling provides 
very efficient ways to control electronic and mechanical degrees of freedom of NEM-devices. 
The conEned area of movable nano-meter sized island of electron gas (quantum dot) is char¬ 
acterized by quantized energy spectrum. The quantum mechanical tunneling between the 
source/drain and quantum dot is responsible for a one-by-one electron charge transfer. Such 
electron transport by periodically moving quantum dot is known as ’shuttling phenomenon’ 
0E]- The signature of shuttling was experimentally demonstrated in Refs. IHT-fTTl 

Recent experimental work [12] suggested a new type of a three-terminal NEM-device as 
a current switch controlled by shifting the frequency of input signal. The ”Y-switch” device 
consisted of three electric terminals and mechanical shuttle component - metallic island 
on top of nano-pillar mechanical resonator H2). The three-terminal device demonstrated 
tunable mechanical modes operating in radio-frequency (RF) regime at room temperatures. 
The applications of this type of NEM-based device for quantum information processing 
potentially include (but not limited to) frequency dependent RF switches and ultra low- 
power logic elements. In our theoretical work we propose an idea of another three-terminal 
device where mechanical resonator plays also a part of one of electric terminals. We suggest 
a mechanism of controlling the switching regime by magnetic field. The high sensitivity of 
NEM resonator provides an opportunity to manipulate the charge transfer in the situation 
when the state of device is defined by out of equilibrium conditions. These systems are 
referred as ’’active NEM devices”. 

In this Letter, we consider a NEM system containing a cantilever as a source located at the 
symmetric point between two vertical drains separated by the air gap of width 2 D as shown 
in Fig. [I] (a). We find that dynamical behaviour of the system shows four distinguishable 
regimes of mechanical vibration as a function of source-drain bias voltage: (i) single stable 
fixed point, (ii) two stable fixed points, (iii) two limiting cycles, and (iv) single limiting 
cycle. Focusing on the regimes (ii) and (iii), we End conditions for the transition between 
two stable states tunable by switching the pulse-shaped magnetic field. 
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FIG. 1: (a) Schematic figure of the system we consider, (b) Equivalent electric circuit (left panel) 
and electromechanical potential U e s(x) in the presence of perpendicular magnetic field (right panel). 
(See the discussion of switching induced by magnetic field in the text). 

We start with formalism of equivalent circuit model for tunnel junction combined with 
Newtonian dynamics of the cantilever of length L vibrating within the air gap width 2D as 
shown in Fig. [I] (b). First, we describe mechanical degrees of freedom of the cantilever by 
the displacement u(z) characterizing the cantilever deflection from the straight configuration 
at the point z (that is, the cantilever axis with its origin at the fixed end, see Fig l.a). 
Second, we introduce the eigenmode representation for the cantilever displacement Pi and 
characterize the fundamental mode of the vibration by the amplitude x. While two position- 
dependent inverse capacitances Cf l {x) of the parallel-plate capacitors are given by a linear 
function of x, the tunnel resistances of the air gap R±(x) exponentially depend on the width 
of the barrier: 

C± (x) = C 0 -p-——, R±(x) = R 0 e* x/X . 

V =f x 

Here, Co (. Ro ) is capacitance (resistance) of air gap of the width D, and A is a phenomeno¬ 
logical tunneling length. The force acting on the cantilever is the vector sum of electrostatic 
force applied to the end of the cantilever, and Lorentz force induced by magnetic field 
B = B 0 e y . The current / through the cantilever and the induced bias across the junction 
V c satisfies the Ohm’s law: / = (V — V c )/R w . We use following notations: R w is resistance 
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of the cantilever, V c = Q c /C(x), Q c is total charge accumulated inside the parallel-plate 
capacitors and C(x) = C + (x) + C_(x). Time evolution of the charge accumulated inside the 
parallel-plate capacitors can be written as: 


Qc 


V 

Qc--b~ = o, 

-Ll'in 


( 1 ) 


Jc(x) T w (x) 

where r c (x) = R(x)C(x ) is a position-dependent RC-time of the tunnel junction, t w = 
R w C(x). Here R{x) = [t/R + (x) + 1/R_(x)}~ 1 = Rq{2 cosh(x/A)) _1 . Then, the electrostatic 
force applied to the cantilever is given by F c = — (Q c ) 2 (dC~ l (x)/dx)/ 2 • e x , and the effective 
Lorentz force induced by the the current is Fb = LI x B. Thus, the equation of motion for 
the vibrating cantilever within the air gap along ^-direction is given by, 


x + 2 , y 0 x + UqX = ai - 


Ql , lb 0 nr IM 

-X + 02-— (V — V c ) , 


( 2 ) 


CoD 2 m ' *7nR 

where m is the effective mass, ojq is a frequency of fundamental mechanical mode and 70 is 
its damping coefficient. Here, 01,02 ~ 1 are geometrical factors. 

In order to present a system of coupled equations describing both mechanical motion and 
charge distribution in compact form, we introduce dimensionless variables, denoted by tilde, 
which are defined by re-scaling the displacement with A, the time with cu^ 1 , the current 
with euj 0 , the voltage with e/C 0l and the force with moj^X [x — x/X, t — u 0 t, q c = Q c /e, 
/ = I/euo, v = CoV/e , and F = F/itiuJqX ): 


~ . 1 ~ . ~ OL ~ 2 ~ , ( ~ ~ 1 

x + —x + x = ~^q c x -\ - (v - q c ) , 

Qo d r 0 r w 

1 v 

q c + - {r w cosh(x) + 1 )q c = —■ 

r 0 r w r 0 r w 


( 3 ) 

( 4 ) 


with r 0 = tuoRoCo, d = D/X, r w = R w /Rq, l/Q 0 = 2 ^ q /ojq. Here 4> B = A LB 0 /(h/e) is 
dimensionless flux through the area of triangle with linear sizes determined by the length of 
the cantilever and amplitude of mechanical vibration measured in the units of flux quantum 
0o = h/e. The dimensionless parameters a, and [3 correspond to the charging energy 
E c = e 2 /C 0 , and oscillator (zero point motion) energy E osc = hco 0 measured in units of elastic 
energy E e = mojQ A 2 : a = a\Ej (mtUgA 2 ), 0 = aoE^J(moj'^X 2 ). Note, that dimensionless 
RC-time To appears in Eqj3] due to rescaling of voltage with the charging energy which is a 
’’natural” unit for rescaling in EqJdJ We ignore the bending effects of the cantilever inside 
the area between source and drain(s), and, assuming that the condition x/d 1 is satisfied, 
disregard the corrections of the order of (x/d) 2 in the equations Eq.([3]) and Eq. (|4j) . 
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Two terms in the r.h.s. of Eq.(|3| account for the retardation effects due to redistribution 
of charge and the Lorentz force acting on the moving cantilever. In the adiabatic limit To 1 
we assume that the dynamics of the charge distribution is determined by RC- time which 
is much faster compared to the time scales responsible for dynamics of mechanical degrees 
of freedom. The approximate analytic solution of Eq. (|4j describes the position-dependent 
charge distribution, q c (x): 


Qc{x) = 


(5) 


r w cosh (re) + 1 

As one can see from Eq. (|5]) , the charge accumulated at the tip of the cantilever decreases 
exponentially with the amplitude x. 

First, let us consider the setup in the absence of perpendicular magnetic fields, B — 0 
(see also [Hj). Then, the effective electromechanical potential, assuming r w <C 1, can be 
written as 




^ x 2 + (In [cosh (x)] 


xtanh(x)). 


( 6 ) 


The values of two local minima x± are found by solving equation d x U e fj(x) = 0. The solution 
reads: x± = icoslW^h-^/a/d 2 — 1 )/r w ). The emergence of x± gives us the condition for 
threshold bias (hi) necessary for formation of a double-well shaped U e g which is written as: 
hi > (r w + 1 )/ \JolI d 2 . 

Nonadiabatic correction to the adiabatic charge 


r = 


hr^ sinh(x) 


• x ■ T 0 


(7) 


( r w cosh(x) + l ) 3 

generates effective ’’negative” friction in the vicinity of the minima of the double-well po¬ 
tential. As a result further increase of bias gives rise to instability of the static state. Finite 
energy pumping generates limiting cycle at the vicinity of two local minima depending on 
the initial condition. With growing bias voltage, the two limiting cycles evolve into one 
large limiting cycle enveloping two local minima. This happens when the energy pumping 
allows the system to overcome the barrier between two local wells. Since there are two char¬ 
acteristic voltages controlling the number of limiting cycles of the system, it is convenient 
to introduce two other critical voltages, V 2 and U 3 for two limiting cycles and one limiting 
cycle, respectively. General expressions for 1 ) 2 , and U 3 are unavailable, however, we evaluate 
the characteristic voltages in Supplementary materials. 
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FIG. 2: (a) Charge distribution q c given by Eq. 5 (black dashed lines) and Eq.([3]) - Eq.Q (green 

dots). The parameter r w = 1, 0.1, 0.01, and, 0.001 (from top to bottom). Other parameters 
are fixed at the values v = 1 and a/d 2 = 1.75. (b) The re-switching current Ijj(black and red 
lines) defined as the current difference between source-drainl and source-drain2 and its fluctuations 
^diif (blue line) averaged over the time interval T = lOOOuyf 1 (two different colors for re-switching 
current correspond to two different initial conditions: black - oscillation near drainl Xi n t = 0.1; red 
- oscillations near drain2, xm = —0.1) as a function of bias. Inset: averaged displacement relative 
to symmetric position of the cantilever and its fluctuations, (c) Poincare map of (x, x) for the 
steady state at zero magnetic field B = 0 evaluated after delay time (t > 5000 <Uq x ) as a function 
of bias voltage. Inset: cross section of main plot at the bias voltages v = 0.8 (black dot), 0.85 
(two red dots) 1.2 (two green curves) 1.6 (blue curve). The parameters are: to = 0.1, /3 = 0.01, 
Q 0 = 100. In order to calculate (a) and (c), we choose 150 random initial conditions in the range 
of inti &inti Qcjint) £ [ 5,5]. 

The main focus of this paper is to demonstrate the re-switching behaviour of active device 
in the situation when the current injected mainly in the drainl is eventually forced to be 
injected to the drain2. This is why we will be interested in the calculation of the difference 
between the current injected from the source to each of the drains, calling the current 
through the drainl as I + and the drain2 as /_ (Fig l.a). The difference between these two 
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currents, which we call a ” re-switching current” I^is, can be fine-tuned by applying external 
magnetic held. The current — /+(t) — J_(f) hilly describes switching properties of the 

active device, ft is convenient to characterize switching by the current averaged over a time 
interval T being large compared to the period of mechanical vibrations (in the numerical 
calculations we use T = 1000 u;^ 1 ). Another important for switching dynamics time scale 
is associated with the delay t after which we perform the time averaging (t > 5000 oyf 1 ) in 
the steady state (see Fig. [2](b) ). For V\ < v < u 3 , the difference I^ s splits following the 
evolution of the system and depending on the initial condition either to the right ( x int > 0, 
black color in Fig. [2](b) ) or to the left ( Xi n t < 0, red color) near one of two minima of the 
potential U e s(x). The shuttling regime (v -2 < v), is characterized by non-zero fluctuations of 
the current difference, — Ifm)- The fluctuation strength continues to increase 

after a sudden drop of I™ s at v 3 , (see the blue line in Fig. |2](b)) . The stationary Poincare 
map of various random initial condition for (x,x) as a function of v is shown in Fig. [2](c). 
The one-to-one correspondence between the re-switching current and displacement shown 
on Fig. [2] (b) can be used for position detection of the nano-device. 

Next, we describe the setup in the presence of perpendicular magnetic held applied in 
order to manipulate the switching current between source and one of two drains, I± in the 
V\ < v < v 3 regime. We consider adiabatically adiabatically varying time-dependent flux 
= ^(tanh(^r^) — tanh(f^)), under following condition for duration of the flux pulse 
Td = t e — t. s and saturation time r p compared to the AC-time r 0 : 1/Qo < T o 'C Td -C r p . 
Fig. 1 (a) to (d) illustrates switching dynamics of a system initially located in the left 
minimum. We apply pulses Td = 200, and t p = 1000 at t s = 15000 {t e = 35000) for 
stimulating jumps from the left well to the right one and back. The voltage dependence of 
the lower/upper critical fields is shown on Fig. [3] (e). If the magnetic held exceeds its upper 
critical limit, the double-well potential transforms into the single-well potential (Fig. [4] 
(b)). The voltage dependence of the upper critical hux cj)g ax can be obtained by evaluating 
the minima of confining potential under condition that two stable minima transfer into 
single stable minimum. The lower critical held have been numerically defined as the 
minimal value of applied hux by comparing x averaged over long time scale T before and 
after stimulation. In Fig. [3] (e), red colored line shows log-scaled current square average, 
(|/_i_(t) + J„(f)| 2 ) at the which is directly proportional to current power. It is therefore 
demonstrated that there exist regimes when small magnetic held can switch large currents in 


7 




V 



FIG. 3: Switching behavior of the device controlled by magnetic pulses. Time scanning of 4>b (a), 
x (b), /+ (c), and /_ (d) at the bias voltage v = 0.85, feo = 1 (black), and v = 1.2, bo = 4 (red), 
(e) Critical magnetic field (black) and (p 1 £ ax (blue) for switching phenomenon as a function 
of bias and the mean-square of the total current averaged at the critical value of magnetic field 
(red), (f) Position-dependent effective friction under different values of dimensionless flux 0 b from 
its minimal valued 0 (black) to maximum value 5 (magenta) with the step A <f>B = 1. 

the active regime of nano-device due to amplification of device sensitivity by the preceding 
signal. As is seen from Fig. [3] (e), small flux switches between two different regimes both 
at the voltages around v\ and v%. This means, that in addition to use of magnetic field for 
manipulating the current switch, one can use switching itself for detection of small magnetic 
fields thus providing a highly sensitive magnetic field sensor. 

The switching mechanism based on magnetic fields in the device can be considered by 
using position-dependent effective dissipation coefficient. Since the device shows position- 
dependent charge distribution, the Lorentz force involves non-adiabatic corrections to charge 
dynamics. Consequently, position-dependent effective dissipation has been emerged in the 
equation of motion: 

'Kj3uJov4>Br w sinh(x) 2au>oTo sinh(x) 

(r w cosh(5) + l) 3 d 2 (r w cosh(x) + l) 4 

Fig. [3] (f) shows position-dependent 7 e ff(x) as a function of varying magnetic field. Unlike 
the Lorentz force, which is an odd function of the coordinate x, flux <pB and voltage v, the 
non-adiabatic contribution to the electrostatic force being even function of both coordinate 
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FIG. 4: Phase diagram for bi-stability in the parameter domain ( v,a/d 2 ) under <f>B = 0 for (a) and 
(cpBi a/d 2 ) with 7 = 1 for (b) by using adiabatic approximation, Eq.([5]). The brown colored area in 
(a) represents the shuttling-promoted switching regime, V 2 < v < v 3 . The details for computing the 
boundaries of each domain are given in the Supplementary materials. Thick lines denote numerical 
solution of Eq. ^ and Eq. Q describing the evolution of critical voltages hi (black), V 2 (red), and v 3 
(blue) as a function of applied bias v. The shape of electromechanical potentials in (b) corresponds 
to the points <j>B = 0 (black), —0.5 (red), and 0.5 (blue) at a = 1.75, and (j>B = —0.75 (magenta), 
0.75 (green) at a = 1.5 

x and bias voltage v always reduces the dissipation near stationary position (see last term 
in Eqj8]). 

Using adiabatic approximation, we calculate the phase diagrams of bi-stability regime (see 
Fig. 0, from which the potentialities of current switch can be seen. We use folowing color 
scheme in Fig. |4|a): gray color is used for ’passive’ switching regime (hi < v < h 2 ), brown 
color denotes the ’active’ switching regime based on the shuttling instability (h 2 < v < v 3 ). 

In conclusion, we have reported current-switching device promoted by shuttling phe¬ 
nomenon based on magnetically controllable bi-stability based on strong NEM coupling. 
The NEM coupling gives rise to double-well shaped electromechanical potential controlled 
by bias voltage between source and drain. Based on the electromechanical pumping provided 
by shuttling phenomenon, the switch can transfer huge current power as an active device. 
We have performed both numerical and analytical analysis and found regime of shuttling 
instability. 
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It worth noting that similar instability and transient from damping NEM oscilations (cool¬ 
ing regime) to self-sustained large-amplitude shuttling (heating regime) may be realized for 
’’Kondo shuttling” [151H5] between metallic lead and long metallic cantilever with attached 
nano-island HZ). It is tempting to unite spin and charge related switching mechanisms in 
the same system, e.g. in a three-terminal device with two Kondo shuttles. 

The principal scheme of device proposed in this letter may be useful not only for switching 
application in electronic circuits but also for testing magnetically fine-tuned multi-stability 
in non-linear system as well as toy model of dynamic critical phenomena in dissipative 
systems. 

We appreciate fruitful discussions with K.-H. Ahn, Hee Chul Park and S. Ludwig. The 
research of KK was partially supported by ISF Grant No. 400/12. The work of RIS and 
LYG was supported in part by Swedish VR. 
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